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Abstract 

Wc consider a reaction-diffusion equation in narrow random channels. We ap- 
proximate the generahzcd solution to this equation by the corresponding one on a 
random graph. By making use of large deviation analysis we study the asymptotic 
wave front propagation. 
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1 Introduction 

In studying the motion of molecular motors we introduced in a solvable model: 
we think of the molecular motors as diffusion particles traveling in a narrow random 
channel. Based on the model suggested in [5], we consider in this paper wave front 
propagation for a reaction-diffusion in narrow random channels. Problems of this type 
naturally appear in the theory of nerve impulse propagation and in combustion theory. 
Our analysis relies on techniques in large deviations similar to that of [31 Chapter 7] 
and [13], [13], [15], [1], [TH Chapter 5]. We shall note that problems of this type are 
mentioned in |H Chapter 7], |8], [7]. It is also interesting to note that similar problems 
are considered in [12j . |16j . |17j but from different points of view. 

Let us first briefly recall the model introduced in [H]. Let h^{x) be a pair of 
piecewise smooth functions with h^{x) — /iq (^) = l-oi^) > 0. Let Dq = {{x,z) : x € 
M, /ig (a;) < z < h^{x)} be a tubular 2-d domain of infinite length, i.e. it goes along the 
whole X-axis. At the discontinuities of h^{x), we connect the pieces of the boundary via 
straight vertical lines. The domain Dq models the "main" channel in which the motor 
is traveling. Let a sequence of "wings" Dj (j > 1) be attached to Dq. These wings are 
attached to Dq at the discontinuities of the functions h^{x). 
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Fig. 1: A model of the molecular motor. 



Consider the union D = Dq |J |J L'j . An example of such a domain D is shown 

Vi=i J 

in Fig.l, in which one can see four "wings" Di,D2, D^, D4. We assume that, after adding 
the "wings", for the domain D, the boundary dD has two smooth pieces: the upper 
boundary and the lower boundary. Let n{x,z) = {ni{x , z) , n2{x , z)) be the inward unit 
normal vector to dD. We make same assumptions as in [5]. 

Assumption 1. The set of points x S M for which there are points {x,z) € dD at 
which the unit normal vector n{x, z) is parallel to the x-axis: n2(x, z) = has no limit 
points in M. Each such point x corresponds to only one point (x, z) € dD for which 
n2{x,z) = 0. 

Assumption 2. For every x the cross-section of the region D at level x, i.e., the 
set of all points belonging to D with the first coordinate equal to x, consists of either 
one or two intervals that are its connected components. That is to say, in the case 
of one interval this interval corresponds to the "main channel" Dq; and in the case of 
two intervals one of them corresponds to the "main channel" Dq and the other one 
corresponds to the wing. The wing will not have additional branching structure. Also, 
for some < /q < ^0 < 00 we have Iq < /i^j'(x) — /iq (x) = lo(x) < Iq. 

Let us take into account randomness of the domain D. Keeping the above assump- 
tions in mind, we can assume that the functions h^{x) and the shape of the wings 
{k = 1, 2, ...) are all random. Thus we can view the shape of D as random. We introduce 
a filtration J^, — oo<s<t<ooas the smallest cr-algebra corresponding to the shape 
of D n {(x, z) : x G We introduce stationarity and mixing assumptions. Let us 

consider some A ^ J-^, —00 < s < t < 00. The set A consists of some shapes of the 
domain D n {{x,z) : x € [s,t]}- Let 6r {r G M) be the operator corresponding to the 
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shift along x-direction: Or (A) € J^lXr consists of the same shapes as those in A but 
correspond to the domain D n {{x, z) : x G [s + r,t + r]}. 

Assumption 3. (stationarity) We have P(^) = P{9r{A)). 

Assumption 4. (mixing) For any A ^ and any B G J'li.r have 



hm \P{AnB) -P{A)F{B)\ 

r— >itoo 







exponentiahy fast. 

For instance, we can assume that there exists some M > such that P{A D B) = 
P{A)P{B) for |r| > M. 

In particular, the mixing assumption implies that the transformation Or is ergodic. 

Here and below the symbols P and E etc. refer to probabilities and expectations 
etc. with respect to the filtration {/'s}_oo<s<i<oo- 

Let D'^ = {{x,£z) : {x,z) € D}. The parameter e > is small. The domain 
models the narrow random channel. Let us consider the following reaction-diffusion 
equation in the domain D^: 



du' _ 1 
'dt ~2 
u^{0, X, z] 



dv 

e _ 



9V 
dx"^ 
= 9{x) 

, 



+ 



dz"^ 



(LI) 



u^{t, X, z) , (t, X, z) G M+ X D'^ . 

Here v is the inward unit normal vector field on dD^] V{x, z) is the velocity field; 
the function /(n^) is smooth and is of KPP type: / G C(°°)([0, 1]), /(O) = /(I), < 
f{u^) < f'{0)u'^ for all G (0, 1), e.g. f{u^) = n^(l — n^). The initial function g{x) > 
is smooth and compactly supported in x G M: ^'(x) G Cq°°^(]R). We notice that the 
initial function g{x) depends only on the variable x G M and is independent of z. 

Alternatively, problem (1.1) can be considered on the domain D with a change of 
variable z i— > z/e. The equivalent problem takes the form 



du' _ 1 
~dr ~ 2 

u^{0, X, z] 



dx"^ 



£ _ 



(1.2) 







dD 



u'^{t, X, z) , {t, X, z) G M+ X D . 



Here u'^ = (j^fiJ^f) — {^'>T'i{^j ^)^'n'2{x, z)) is the inward unit co-normal vector field 

1/92 1 92 

on dD corresponding to the operator - -— r H — ^tt^ 

2 \ax^ az^ 

The diffusion process Xf = {Xf , Zf) corresponding to problem (1.2) takes the form 
of equation (2) in [5]. We have 
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dx! = dwl + v{xt, zt)dt + viixf, zi)di^ 



t ) 



e 

Let P'^, denote probabilities and expectations with respect to the filtration 
generated by (and henceforth {Wl,Wi)). We note that as in [5] the motion of 
is independent of the random shape of D (and henceforth D^). We have, in the same 
way as in [5], the following. 

Assumption 5. The process iW^, Wf) is independent of the filtration {Fl} ~oo<s<t<oo 
corresponding to the shape of D. 

We shall also make some assumptions parallel to Assumptions 7 and 8 in [5]. To 
this end we let L be the random variable distributed the same as the distance along 
X-axis between two wings: L is the distance along x-axis between two cross-sections of 
D where there is a branching. Let Zwing be the cross-section width of the wing. Let r 
be the projection length of a wing onto x-axis (r can be positive or negative; compare 
with O Assumptions 7 and 8]). We assume the following. 

Assumption 6. With P probability 1 we have (1) -|-oo > L > L > L> for some 
constants L,L > 0; (2) < /wing(a;) < Ai and \r\ < Ai for some constant Ai > 0. 

Our goal in this paper is to study the asymptotic wave front propagation properties 
for the generalized solution of (1.2). To be precise, by a generalized solution of (1.2) we 
mean the one defined via the path integral representation (Feynmann-Kac) formula: 



u''(t,{x,z)) = E 



(x,z) 



exp (^l\{u'{t-s,Xl))ds^ g{Xt 



(1.3) 



Here c(n) = for u > and c(0) = lim = sup ^ . The latter equality is 

due to the KPP nonlinearity assumption. We shall also suppose that |c'(u)| < Lip(c) < 
oo, u € [0,1]. The proof of existence, uniqueness and regularity of the generalized 
solution to the integral equation (1.3) is close to [3l Chapter 5, Section 3]. For the 
reader's convenience we will prove it in Section 3 of this paper. 

We introduced in [5] the metric graph F corresponding to the domain D (see Fig.l). 
Let Yt = {Xt,kt) be defined on F as in [5l Section 2]. The construction of the graph 
F and the process Yt as well as some basic convergence results in [5] will be recalled 
in Section 2. The process 1^ is a diffusion process on F with a generator A and the 
domain of definition D(A). We consider the reaction-diffusion equation associated with 
the Markov process Yt. This equation takes the form 



— = Au + f{u) , n(0, {x, k)) = g{x) , n = u{t, y), {t, y) € M+ x F . (1.4) 
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The initial function g = g[x) is the same as in (1.2). We require that for each fixed 
t > 0, u(t, •) € D{A). This requirement ia a kind of boundary condition. For details we 
refer to [5] and [21 Chapter 8] and the references therein. 

The generalized solution to (1.4) is defined as the solution to the integral equation 



u{t,{x,k)) = E(^^fc) 



exp c{u{t - s,Ys))ds^ g{Xt 



(1.5) 



The existence and regularity of the solution can be proved in a same way as those 
for (1.3). We will briefly mention this in Section 3. 

We will show, in Section 3 of this paper, that as e is small, the solution u^(t, {x, z)) 
of (1.3) will converge to u{t,{x,k)) of (1.5) in the strong sense. Here {x,k) = 2)((x,z)) 
and the mapping 2) : D — )• T is an identification map that will be recalled in Section 2. 

After we get convergence results we will focus on the study of the solution u{t, (x, k)) 
of (1.5). We will show that, as f — > oo, the solution u{t, {x,k)) behaves asymptotically 
as a traveling wave. This wave is traveling in both positive and negative directions along 
X-axis. 

The paper is organized as follows: in Section 2 we recall some necessary basic set 
up of [5], such as the construction of the graph T, the process It, etc.; in Section 3 we 
show the convergence as e is small of the solution u^{t, (x, z)) to u{t, (x, k)); in Section 
4 we obtain some auxiliary results that will be used in later sections; in Section 5 we 
derive the large deviation principle; in Section 6 we study the wave front propagation 
properties of the solution u{t, (x, k)). 



2 Set up 

We shall first recall some basic facts in [5] (also see [TO]). Let us work with a fixed 
shape of D. 

First of all we need to construct a graph P related to the domain D (see Fig.l). 
For xo € M, let C{xq) = {(x, z) ^ D : x = xq} be the intersection of the domain D with 
the line {x = xq}. The set C(xo) may have several connected components. We identify 
all points in each connected component and the set thus obtained, equipped with the 
natural topology, is homeomorphic to a graph P. We label the edges of this graph P by 
/i, /fc, ... (there might be infinitely many such edges). 

We see that the structure of the graph P consists of many edges (such as Ii, I3, 15, Ij, Ig,... 
in Fig.l) that form a long line corresponding to the domain Dq and many other short 
edges (such as I2, h, h, Is,--- in Fig.l) attached to the long line in a random way. We 
will henceforth denote by Iq the long line corresponding to the domain Dq. 

A point y G P can be characterized by two coordinates: the horizontal coordinate 
X, and the discrete coordinate k being the number of the edge in the graph P to which 
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the point y belongs. Let the identification mapping be 2) : -D ^ T. We note that the 
second coordinate is not chosen in a unique way: for y being an interior vertex 0/ of 
the graph F we can take k to be the number of any of the several edges meeting at the 
vertex Oj. 

The distance p{yi,y2) between two points yi = {xi,k) and y2 = {x2,k) belonging 
to the same edge of the graph T is defined as p{yi-,y2) = l^^i — a^2|; for yi,y2 G T 
belonging to different edges of the graph it is defined as the geodesic distance p{yi,y2) = 
min(/?(yi , Oj-^ ) + p{Oj-^ , Oj^ ) + ... + p{Oj^ , 2/2 ) ) i where the minimum is taken over all chains 
yi -f-)- Oj^ o -f-)- ... o Oj, -H- y2 of vertices Oj. connecting the points yi and y2- 

For an edge Ik = {{x,k) : < x < B^} we consider the "tube" Uk = 2)~^(/fc) H 
{Ak ^ X < Bk] in D. The "tube" can be characterized by the interval x G [^^,5^] 
and the "height functions" h^{x): = {{x,z) ■ A^ < x < Bk,h'^{x) < z < K^{x)}. 
For X € [Ak-, -Bfe], we denote the set Ck{x) to be the connected component of C{x) that 
corresponds to the "tube" Uk'- Ck{x) = {x} x Let lk{x) = h'^{x) — h~^{x) > 

for all X G R. We notice that each h^{x),lk{x), etc. is smooth. 

The vertices Oj correspond to the connected components containing points (x, z) G 
dD with n2{x^z) = 0. There are two types of vertices: the interior vertices (in Fig.l 
they are 01,03,05,07) are the intersection of three edges; the exterior vertices (in Fig.l 
they are O2, O4, Oe, 0%) are the endpoints of only one edge. 

Using the ideas in [lOj with a little modification we can establish the weak conver- 
gence of the process = 2)(JCf ) (which is not Markov in general) as e 4, in the space 
C[o,r](r) to a certain Markov process on F. A sketch of the proof of this fact is in [5l 
Section 2]. 

The process Yj is a diffusion process on F with a generator A and the domain of 
definition D{A). We are going now to define the operator A and its domain of definition 
D{A). 

For each edge Ik we define an operator L^: 

— 1 d f du\ — du 

^'^^("^ = ^)d^ r J + ^^-("^ ^A>^^-^Bk. 

Here 

_ 1 /■'^feW 

^k{x) = Y-r^j V{x,z)dz 

'■k{X) Jh-{x) 

is the average of the velocity field V{x,z) on the connected component Ck{x), with 
respect to Lebesgue measure in z-direction. At places where Ik = 0, the above expression 
for Vkix) is understood limit as Ik — )• 0: 

_ 1 fht(y) 

Vk{x) = hm — - / V{y,z)dz . 
lk{y) Jh-iy) 

We will assume throughout this paper the following. 
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Assumption 7. The function Vk{x) = 0. 

We notice that this is a bit different from the corresponding one in [5l Assumption 
6]. We point out that the vanishing mean drift assumption is crucial for the method of 
our analysis to work. 

Thus under our Assumption 7 we have 

The operator Lk can be represented as a generalized second order differential op- 
erator (see [2]) 

Lku{x) = Dm.^DpJ{x) , 

q{x -\- 51 — (j(x^ 

where, for an increasing function h, the derivative is defined by Dhg(x) = lim — — , 

610 h[x + 0) — h{x) 

and 



Pk{x) 



dx 



h{x) 

n^-kix) = 2 / lk{x)dx 



is the scale function, 



is the speed measure. 

The operator A is acting on functions / on the graph T: for y = (x, k) being an 
interior point of the edge we take Af{y) = Li^f{x, k). 

The domain of definition D{A) of the operator A consists of such functions / 
satisfying the following properties. 

• The function / is a continuous function that is twice continuously differentiable 
in X in the interior part of every edge 1^', 

• There exist finite limits lim Af{y) (which are taken as the value of the function 

Af at the point Oj); 

• There exist finite one-sided limits lim Dp f{x,k) along every edge ending at 
Oi = {xi, k) and they satisfy the gluing conditions 

5^(±) liml?p /(x,%) = 0, (2.1) 

^ * X^Xi J 

where the sign "+" is taken if the values of x for points (x, kj) € Ikj are > Xi and "— " 
otherwise. Here Ni = 1 (when Oi is an exterior vertex) or 3 (when Oi is an interior 
vertex) . 

For an exterior vertex Oi = (xi, k) with only one edge Ik attached to it the condition 
(2.1) is just lim Dp^f{x,k) = 0. Such a boundary condition can also be expressed in 
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d ^ f 

terms of the usual derivatives — instead of Dj,, . It is lim lk(x)-—(x,k) = 0. We 

dx x^xi ^ ' dx 

remark that we are in dimension 2 so that these exterior vertices are accessible, and the 
boundary condition can be understood as a kind of (not very standard) instantaneous 
reflection. In dimension 3 or higher these endpoints do not need a boundary condition, 
they are just inaccessible. For an interior vertex the gluing condition (2.1) can be 

written with the derivatives — instead of Dp, . For k being one of the kj we define 

dx 

= lim lk{x) (for each edge 1^ the limit is a one-sided one). Then the condition (2.1) 
can be written as 

^(±)a,,^, . lim = . (2.2) 

^-^ x^Xi dx 

i=i 

It can be shown as in jlO| Section 2] that the process Yt exists as a continuous 
strong Markov process on F. 

We fix the shape of D. For every e > 0, every x = {x,z) ^ D and every T G (0, oo) 
let us consider the distribution of the trajectory = 2)(JCf) starting from a point 
Xq = ic in the space C[o_r](F) of continuous functions on the interval [0,T] with values 
in F: the probability measure defined for every Borel subset B C C[o^r](F) as fJ.%{B) = 
P^£:^^(y/ G B). Similarly, for every y G F and T > let fiy be the distribution of the 
process Yt in the same space: fiy{B) = P^(y, G B). The following theorem is basic for 
our analysis. 

Theorem 2.1. For every x ^ D and every T > the distribution fj,^ converges 
weakly to to Aig(a.) as e J, 0. 

In other words we have 



for every bounded continuous functional F on the space C[o,r](r)- 

The proof of this theorem follows from [TO] and there is a sketch in [SI Section 3] . 
We omit duplicating the details here. 

3 Convergence of to u 

We recall that our definition of the generalized solutions to (1.2) and (1.4) are the 
solutions of the integral equations (1.3) and (1.5), respectively. That is to say, we have 



n^(t, {x,z)) = E 



ix,z) 



exp(^l\{u%t-s,X2)ds^ 9{X!) 
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and 



(x, A:)) = E(3.^fc) exp^y c{u{t - s,Ys))ds^ g{Xt) 



Theorem 3.1. There exist unique bounded measurable generalized solutions u'^{t, {x, z)), 
t > 0, [x, z) £ D and u{t, {x, k)), t > 0, (x, /c) € F for (1.2) and (1-4), respectively. These 
solutions are continuous for all t > 0. 



Proof. We take (1.2) as an example. The proof for (1.4) is exactly the same. 
We shall prove the existence and regularity by using a contraction mapping principle 
(compare with [3l §5.3]). To this end we consider the Banach space Bj- of bounded 
measurable functions on [0,T] x [D] with norm ||f || = sup |t;(t,x)|. Consider in 



Bt the following operator 



{t,x)€lO,T]x[D] 



F[v] = F[v]it,ix,z))=B^^^,) 
It is then checked that we have, for < t <Tq <T, that 



exp (^l\{v{t-s,X',))ds^ g{Xl) 



V e Bt 



u]-F[v]\{t,{x,z)) 



E 



{x,z) 



E 



< 



E 



{x,z) 



{x,z) 

t 



exp c{u{t-s,Xl))ds^ g{Xl) 

exp c{u{t — s, Xl))ds-^ — ex.p c{v{t — s, Xl))ds 

< \\g\\ exp(c(0)t) • Lip(c) ■ t ■ \\u — v\\ 
= C{To)\\u-v\\ 



exp (^l\{v{t-s,Xl)ds^ g{Xt) 



for u,v € Bt and C{Tq) 



contraction provided that Tq < 



\\g\\ exp(c(0)ro)Lip(c)ro. 
1 



This guarantees that F is a 
By contraction mapping the- 



|c/|lLip(c) exp(c(0)r) 

orem we have existence and uniqueness of generalized solution in the space of bounded 
measurable functions to the problem (1.2) on the interval [0,To]. Since the solution 
\u^{Tq, (x, z))\ < exp(c(0)To) we can use the same Tq and work with intervals 
[To,2To], up to [{k — l)TQ,kTo] {k > 1) provided that kTo < T. This gives exis- 
tence "in the large" for a unique generalized solution u^{t, {x,z)) for (1.2) in the space 
of bounded measurable functions. 

The continuity of the solution u'^{t, (x, z)) in the variables t and (x, z) is provided 
by (1) The continuity of g; (2) The Lipschitz continuity of c(u); (3) The continuity and 
continuous dependence of the process Xf on Xq = (x, z). □ 
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We shall then show the approximation of the generalized solution u^(t,{x,z)) as 
e is small by the generalized solution n(t, 2)((x, z))). We prove this via a sequence of 
auxiliary results. 



Lemma 3.1. We have 



lim max max \u^(t, (x, zi)) — u'^it, (x, Z2))\ = ■ (3-1) 

£lO (x,2i),(x,22)e2)-i((a;,fc)) 0<t<T 



Proof. We consider the stopping time = {{x , Z2) , zi) = inf{t > : X; 
X, Zq = Z2, Zf = zi\. By strong Markov property of the process we have 



E 



E 



{x,Z2) 



exp ^ c{u'{t - s, XI 



exp( / c{u%t-s,Xl,Zl))ds]u%t-T',{X^^e,zi)) 



Since the motion is moving very fast as e J, we see that ^ almost surely 
as e 4, 0. This immediately implies the convergence. □ 



Let {x, k) € P. We introduce a new function 



u'it,ix,k)) 



|2) ^{{x,k))\J<Q^i(^^^^k)) 
and we see from Lemma 3.1 that we have the following 



u^{t, (x, z))dz 



Corollary 3.1. We have 



lim max max \u'^(t,^)(x, z)) — u^(t, (x, z))\ =0 . 



(3.2) 



We are going now to prove that the function u^{t, (x, k)) has a uniform in e bounded 
first derivative in the variable x. 



Lemma 3.2. We have an a-priori estimate 



max 

o<t<T,{x,k)er 



dx 



(t, {x,k)) 



< C(T) = C 



(3.3) 



where C > is independent of e. 



Proof. By (1.3) we have 
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u%t, {x,z)) = E(^^^) 



Differentiating with respect to x we have 



exp (^|\{u'{t-s,Xl))ds^ g{X.^ 



dx 



+<7(Xf)exp(^^ c{u'(t-s,XlZt))^-^(^J^ ciu%t-s,XlZl))ds 



Note that 



dx \./o 



c(u^(t - s,Xl,Zl))ds 



c'iu^t - s,Xt, Zl)) — {t - s,X!,Zt)^ds . 







dx 



Therefore if we let rrfit) = max 

{x,k)(^D,0<s<t 



^{s,{x,k)) 



rrf{t) < a{t) + /3(t) / rrf{s)ds , 



, we get 



where a{t),j3{t) are bounded with their bound depending on the regularity of g{x), 
c{u) and the shape parameter /(x), yet independent of e. We then apply a Gronwall 
inequality to conclude that we have an a-priori estimate 



max 

0<t<T,{x,z)&D 



dx 



{t, {x,z)) 



< C{T) = C 



(3.4) 



where C > is independent of e. 

Prom the above estimate and taking into account the smoothness of the shape 
parameter /(x), we see that the a-priori estimate in the statement of the Lemma holds. 
In fact, by the definition of u^(t, (x, /c)), we have 

\^-\{xM\-n'{t,{x,k))= I u'it,ix,z))dz . 



Thus 



// (fll^ ft 

— {\^-\{x,k))\W{t,ix,k))+W-Hix,k))\ — {t,{x,k)) = — 



dx 



dx 



2]-i({a;,fc)) 



u^{t, (x, z))dz 



We notice that 



d 



u^{t, (x, z))dz 



dx yj<g-l({a;,fc)) 

Here ^~^{x,k) = (a(x),6(x)). This implies (3.3). □ 



du^ 

— — (t, (x, z))dz+u'^{t, (x, b{x)))b'{x)—u'^{t, (x, a(x)))a'(x) 

^-'^({x,k)) dx 
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Making use of Theorem 2.1, Corollary 3.1 and Lemma 3.2 we can prove the follow- 
ing. 



Theorem 3.2. We have 



lim max max \u^{t, (x, z)) — u{t,^X){{x, z)))\ = . 
elO 0<t<T (x,z)eD 



(3.5) 



Proof. An outline of this proof is mentioned at the end of [U Chapter 7]. We 
fulfill the details here. 

Let {x,k) = 2)((x,2;)). From (1.3) and (1.5) we see that 

|u^(t, {x,z)) - u{t, {x,k))\ 



E 



{x,z) 



exp c{u%t - s, Xl))ds^ g{Xl) - E(^^fc) exp (^j^ c{u{t - s, Ys))ds^ g{Xt 



< \\g\\ ■ [(/) + (//) + (///)] . 
Here 



(I) 



E(.,.)exp(^^ c{u^t-s,X'^))ds^ -E(,,,)expQ c{u%t - s,^{Xl))ds 



{II) 



E(x,,)exp(^^ c{if{t-sMXl)))d^ - E(,,fc) exp (^^ c{u' {t - s,Ys)d^ 



{III) 



E 



x,fc)exp c{u^t - s,Ys))ds^ -^(^^^k^exp (^j c{u{t - s,Ys)ds 



Thus we see that (!) ^ as e J, due to (3.2) ; (//) — )■ as e J, due to the weak con- 
vergence of the processes Yf = 2)(Xf ) to It on P as e i in C[o,t] (r) (Theorem 2.1) and 

(3.3); (///) can be bounded by a constant multiple of / max \u^{s,{x,z)) — u{s,{x,k))\ds 

Jo (x,z)<^D 

plus a term going to as e J, (due to (3.2)). We then apply a standard technique via 
Gronwall's inequality and we can conclude. □ 

4 Auxiliary results 



This section will be devoted to obtaining some auxiliary results which will be used 
in Sections 5-6. 
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Let the random variable 



= inf{t >0,Yo = (s, 0),Yt G Jq, < r} (4.1.1) 

for s > r € M. Intuitively, is the first time that the process Yt, starting from 
Yq = (s,0), comes back to Iq with the value of its x-component < r. We recall that Iq 
is the long line in T corresponding to the domain Dq. 
In the same way we define 

= inf{t >0,Yo = (s, 0), Yt G Iq, Xt > r} (4.1.2) 

for s < r G M. 

Let A > 0. Let the function 

n(x) = E^[e-^^o] (4.2) 

for x G R. We remind the reader of a small notational convention here. In this section 
for convenience of notation we have a minus sign in front of the stopping time Tq in 
(4.2). In the Sections 5-6 we will drop this minus sign and instead we will be mainly 
working with A < 0. 

Let us first consider the case when x > 0. The function u{x) is the solution of the 
following Sturm-Liouville problem 



- Au = on r , u G D{A) , m(0) = 1 , m(+oo) = . (4.3) 

To solve the above problem we shall first recall the basic theory of Feller ([2]). We 
follow here [11] and we also refer the reader to [6l Lemma 2.10]. Let first some r > 0. We 
consider the eigenvalue problem associated with the generalized second-order differential 
operator D„iDp: 

DmDpu{x) - Xu{x) = (4.4) 

on an interval x G / = [0, r]. Here m = m(x) is the speed measure and p = p{x) is the 
scale function. The function p(x) is a strictly increasing continuous function on (0, r) 
and the function m(x) is a strictly increasing function on (0, r) continuous to the right. 
The generalized derivatives are defined as 

Dpf{x) = hm — -— , Dmf{x) = hm 



y^x p[y) — p[x) ' y^x m{y) — m{x) 

where x G (0, r) and / is a real function defined in a neighborhood of x. There are 
two basic solutions u+(x), u-{x) of the equation (4.4) with tt+(0) = U-{r) = and 
U-)_(r) = U-{0) = 1; the function u^{x) is increasing in x and ii_(x) is decreasing in x; 
the derivatives DpU^{x), DpU^{x) are increasing functions. 
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Moreover, an explicit representation of the functions u±{x) is available ([llj). 
set, for n = 0, 1, 2, x G [0, r], 

u^'^^x) = 1 , u'^^+^Xx) = [\^''\s)dm{s)dp{y) . 

Let 



JO 



(a;,A) = ^A"u(")(x) . (4.5 



n=0 



It could be justified that the above series converges. 
We can easily check that 



n=0 

Let 



Dpu{x, A) = V A"+i r u(") iy)dmiy) . (4.5. 



u+{x,X) = u{x,X) / u{y,X) "^dpiy) , 



U-{x,X) = u{x,X) / u{y,X) '^dp{y) . 

J X 

Moreover, we can calculate the derivative 



dx ^ dx ^ Jq ' u{x, A) dx 



' dx ^ J ^ ' m(x, A) 

Then we have 

n+(x,A) 'u_(x,A) 

u+(x) = ; — — , n_(x) = -. — -T, (4.5 

+ ^ ' ii+(r,A) ^ ' u_(0,A) ^ 



and 



duj^ , , _ dx ' 7o ' ''^(^j A) dx 

^ u(r, A) / u{y,X)~'^dp{y) 

r uiyAr^dpiy) ' 'P^^^ 



du- ^ ^ dx ^ ' u{x, A) dx 

n(0,A) r u{y,X)-^dp{y) 



In terms of generalized derivatives we see that the above is equivalent to 



(4.5 



14 



Dpu{x,X) I u{y,X) '^dp{y) + 







u{x, A) 



DpU+{x) = 

u{r,X) I u{y,X)-'^dp{y) 

1 (4-5-2') 

Dpu{x, A) / u{y, A) dp{y) — — 

Dpu^ix) = L . 

u(0,A) / u{y,X)-^dp{y) 
Jo 

A general solution of (4.4) can be represented as a linear combination 

u{x) = C^li+(x) + C~U-{x) . 

The constants and c~ are determined by boundary conditions to be specified. 

In the case when r < situation is similar and we have to make small changes 
accordingly. To be more precise, we can treat the point r as the point and the point 
as the point r. The formulas (4.5.0) ((4.5.0')), (4.5.1) and (4.5.2) ((4.5.2')) have to be 
changed accordingly. In the rest of this section we will be mainly performing detailed 
steps in the calculations assuming r > and we will present corresponding results when 
r < without a detailed calculation. 

Let us come back to our problem (4.3). First of all we note that the structure of 
the graph T consists of two types of edges: the first type of edges are lined up together 
forming the edge Iq and we label them as hk-i, A; € N; the second type of edges 
correspond to the wings and we label them as /2fc, /c S N. These edges are labeled in 
a consecutive way (see Fig.l). Let the projection of the second type of edges l2k onto 
the x-direction be isomorphic to [0,r2fc] for r2k > and [^2^,0] for r2k < 0. Let the 
interval l2k~i be isomorphic to [0, r2k-i]- We solve the problem L^Uk — Xuk = on each 
edge hk-i = [0,r2fc_i] (the first type), hk = [0,r2fc] for r2k > and hk+i = [r2k,0] for 
'''2k < (the second type). We notice that in this case when we represent the operator 
Lfc as a generalized second order derivative operator = Dm^Dp^ we will have 

dmk{x) = 2lk{x)dx (4.5.3) 

and 

dpk{x) = -r-r-^dx . (4.5.4) 

The general solution is represented as Uk{x) = c^u+(x) + c^n^(x). Here u\{x) 
and u^{x) are the two basic solutions corresponding to the interval Ik and we identify 
X with some x £ (or its projection onto the x-axis, anyway). We note that they are 
random solutions. The constants and are to be determined. We shall seek for a 
solution u{x) = U2k-i{x) whenever (x,0) G T. Thus we have 
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U2k-l{r2k-l) = U2k+l{0) = U2k{0) , 

dU2k-l , \ ■ I \ d'^'i.k fr,\ , a dU2k+l .„\ 
a2k-i — j—{r2k-i) = sign(r2fc)72fc-i— (0) +/^2fc+i— ^ (0) 



(4.6) 



OX ax ax 

^xi(O) = 1 , 

lim n2fc+i(r2fc+i) = . 

, fc— )-oo 

In the above a2k~i, f32k+i,72k are the corresponding cross-section width of the 
channel at the junctions. We have a2k-i — /32fc+i = sign(r2fc)72fc. 

Lemma 4.1. We have 



'^2k-l 



-2k~lj 



(4.7) 



Here 



with 



' 1 ' 

, Xk Vk , 



Xk 



hk+l 



dx 



(0) 



dx 



-{r2k- 



(4.8) 



du 



2k -1 



Vk 



dx 



-{r2k~ 



du 



2k~l 



dx 



-{r2k-i) 



+ 



l2k 

a2k-i 



7 dul\ , 



Dp^kKi'^2k) du 



2k-l 



du 



2k 



+ 



dx 



-(0) 



' ,2fe+l 



dx 



-{r2k~i) 



du 



2k~l 



dx 



-{r2k^ 



du'' 



+ 



/32fc+i dx 



-(0) 



a2k-i du''' 



2k-l 



dx 



-{r2k-i) 



(4.9.1) 



if r2k > 0; and 
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du 



2k-l 



Vk 



dx 



-{r2k-i) 



du 



2k-l 



-{r2k-i) 



l2k 
"2fc-l 



dx, 
( dn^ 

dx 



(0) 



du 



2fc-l 



du 



¥+1 



hk+i dx 



■[r2k-ij 
(0) 



"2A;-1 du 



2k-l 



-ir2k-l) 



dx 



{r2k-i) 



dx 

ifr2k < 0. 

Proof. The first three equahties in (4.6) will give us 



(4.9.2) 



"2^-1 
^2fc-l> 



-^2k+l 
^2k+l/ 



where we can calculate the random matrix M^,. Let us first consider the case when 
^2fc > 0. From the first equation of (4.6) we see that we have 



4kDp,,ufir2k) + C2kDp2,U-Hr2k) = . 



(4.10) 



We have ^i2fc- 1(^2^-1) = c^k~i^f H^2fc-i) = <4k-v U2k+i{0) = C2fc+iU^''"^^(0) 



2fc+l/ 



'^2fc+i' ^2fc(0) = C2^u_ (0) = C2^. So from the second equality of (4.6) we get 



'^2k-l ~ '^2k+l ~ ^2k ■ 



(4.11) 



The third equality in (4.6) gives us 



L du+ \ . - du 

a2k-l 1 C^k-l (^2fc-l) +C2fc„i^ 



2fc-l 



{r2k-ij 



du^_t 



du 



2k 



du^_t+^ 



72fc 1 4fc^(0) + C2"fc^(0) j + P2k+i +i^^(0) + 
Prom (4.10), (4.11) and (4.12) we can conclude that 



du 



2k+l 



dx 



-(0) 



(4.12) 



\Xk Vk J 

Here the random variables x^ and are defined by (4.8) and (4.9.1). 
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In the case when < we just have to replace in (4.11) the coefficient by c^^ 
and we have to change the sign in front of 72^ in (4.12) to minus. We thus get (4.9.2). 

□ 

k 

Let Lfc = ^ ^2j-i- We see that without loss of generality we can assume that Li 
and L have the same distribution. 

Lemma 4.2. The random variables — — , /c G N are identically distributed. 
Proof. We have, by strong Markov property of the process Yt on T, that for A > 0, 

By stationarity we see that E^[e ^fe ] has the same distribution as E^[e~'^"^o ]. 
Since 

E^[e-^^<^'^"^]=n(L.-,i) = 4_, 

and 

E^[e-<'=] = n(Lfc) = c-,_i, 
c+ 

we see that the distribution of — — is independent of /c € N. □ 

Lemma 4.3. For any k G N we have < '^'^^ < i almost surely with respect to 



Proof. This is because we have E'^ie-^^o 1 1 e [0, 11 and Lemma 4.2. □ 



It is convenient to introduce the notation 

S'fc = sign(rfc) / Ukiy.XY'^ '^^ 



Here Ufc(x, A) is given by the formula (4.5.0) specified in the interval Ik- 
Lemma 4.4. For any k we have 

_ U2k-l{r2k~l, A)52fc-1 
U2k+l{r2k+l, X)S2k+l 

yk > almost surely with respect to P. 
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Proof. By Lemma 4.1 we have 



2k-l 



'^2k-l ^2k+l , 

Xk— \- Vk 

^2k+l 

Making use of (4.5.2) and (4.5.3) we see that Xk = i ( 2fc 1 1 ) 2k i ^ 'pj^^^g 

U2k+l\r2k+l, A)62fc+1 

rcfc < and since we have Lemma 4.2 and Lemma 4.3 we see that yt > 0, almost surely 
with respect to P. □ 

Lemma 4.5. We have — < < 1 and thus > 1 almost surely with respect 

Vk 

to P. 



Proof. This is a simple consequence of Lemma 4.1, Lemma 4.3 and Lemma 4.4. 

□ 

Lemma 4.6. We have 

1 dU2k-l, 

U2k~i{r2k~iA) dx 

, >N 5*2^-1 U2k{_r2k,^)Dp^,^U2kir2kA)S2k , / ^x5'2fc_l 

+-"2^-1(^2^-1, A)— — — — + ■U2fe-l(r2fe_l, A)- , 

^2k U2k{r2k, AjiJp2fc^2fc(?'2fc, A)62fc + 1 02fe+l 

(4.13.1) 

when r2k > 0; and 



1 , dU2k~l( 

7 rr + «2A;-1 J ('^2fe-l, A)62fe-1 

^^2fc-l('^2fc-l, A) dx 

M2fc-l(r-2fc-l, A) du2fc „ , , ..S2k-1 
+l2k JTTW 1— (0' A)S2fc-l + U2fc_l(r2fc_l, A)- 

U2fc(0,A) dx b2k+i 

when r2k < 0. 



(4.13.2) 



Proof. Let us first consider the case when r2k > 0. By (4.9.1) we can write 



l2k , xs , hk+ 

Vk = yt) -r - 

Here 



Vk = {i) + ^^{{n){iii) + {iv)) + l^{v) . 

02k-l OL2k-l 



dul'~\ ' D,,,uf{r2u) ' dul ^ ^ 

^(0) ^^^(0) 



du^ dut , 
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By making use of (4.5.2) and (4.5.2') as well as (4.5.0') it is straightforward to 
calculate that 



/.N, dU2k-l I ..r, I 

(l) = a2k-l ^ (?'2fc-l, A)52fc-1 + 



{ii) = 
{in) 



ax tt2fc_i(r2fc-i, A) 

U2k{r2kA) 

U2kir2k, A)-Dp2fc'"2fc('^2fc, A)S2fc + 1 

a2fc-l U2k-l{r2k-l, A)5'2fc-1 



72fc U2k{r2kA)S: 



2k 



N OL2k-\ I >N'52fc-l 

{%V) = U2k-l{r2k-lA)—^ , 

72k J 2k 

I N OL2k-\ , ,.S2k~l 

[v) = -o M2fc-i(r2fc_i, A)-^ . 

P2fc+1 '-'2fc+l 



So we get 



"2fc-l — 1 \T2k~\, A)^)2fe-1 + 



dx ' lf2fc-l(?'2fc-l, A) 

1 / .^S2k~\ 
U2k-l{r2k-lA)—F, 



W2fc(r2fc, A)-Dp2fc^i2fc(?'2fc, A)5'2A; + 1 ' S2k 

, . >^5'2fc-l . ,.S2k-l 

+U2k-l{r2k-l, ^ U2k~l{r2k-1, ^)t, • 

J2k J2k+1 

Thus we get (4.13.1). The equality (4.13.2) is obtained in a similar way. □ 

Making use of Lemma 4.6 and basic calculations (4.5.0)-(4.5.4), as well as our 
Assumptions 2 and 6, we see that we have the following. 

Corollary 4.1. 

E In yfc < oo . 



Combining Lemmas 4.5 and 4.6 we colculde that E| lnc^| < Elny^ < oo. Thus we 
have the following. 

Theorem 4.1. 

E[|lnE^[e-^^'f]|] < oo ; E[| lnE^[e-^^o"^]|] < oo . 

In the above theorem the second inequality is estimated in a similar fashion as the 
first one. 
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Lemma 4.7. We have (Tq < oo) = 1 and P^(Tq ^ < oo) = 1 almost surely 
with respect to P. 

Proof. We take Tq as an example. The case for Tq^ is similar. Let v{{x, k); A) = 
P^(tIq'^^ < oo) where {x,k) G F. Here T^q'^^ is the first time the process Yj, starting 
from Yq = {x,k) € T, hits the point (0, /c = 0) G T or {A,k = 0) G P. We set 
Vk{x) = v{{x,k); A) where we identify {x,k) G T with some x G 1^ (or projection of 
onto the x-axis, anyway). Then we have 

lim l2k{x)'^^{x) = , 

dv2k~i, s ■ t \ dv2k,^. . o dv2k+l,^s 
«2fc-i {r2k-i) = sign(r2fc)72fc-^(0) + /j2fc+i (0) , 

■y2fc-l(?^2fc-l) = V2k+l{^) = V2fc(0) , 

t;i(0) = 1 , 
i;((AO);A) = . 

This gives 



V2k{x) = W2fc(0) = ?^2fc-l('^2fc-l) = ^^2fc+l(0) , 



-{x) 



dx hk-lix) 

C2fc+1 = C2k-1 ■ 



Thus 



dy 

V2k+l{x) = D2k+1 + C — - 

Jo Hy) 

We see that Di = 1. And we have recursively that 

r-2fe-l 







Thus we see that 

/o ^o(y) 



vi{x,0);A) = 1 



dy 



A 







dy 
lo{y) 



Thus lim v{{x,0); A) = 1 by Assumption 1 {oo > /q > ^0(2;) > > 0). In particular, 

A— >oo 

pW ^rpL ^ i^-j _ almost surely with respect to P. □ 
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Lemma 4.8. We have 'E^[Tq] = +00 and E^[Tg ^] = +00 almost surely with 
respect to P. 

Proof. We take Tq as an example. The proof for T^^ is similar. We show by 
comparison. To this end we construct the part of the process Xf = {Xf, Z^) within the 

domain Dq (compare with O Section 4.2]). Let ipt = I € DQ)ds be an additive 

Jo 

functional, which is called the proper time of the domain Dq. We introduce the time j3t 
inverse to ipt and continuous on the right. Let Y^'^° = 2)(JC^J. One can show that 
as e 4 the weak convergence of y/'^" to Y^^\ The process is described as a 
one-dimensional diffusion process on M with gluing conditions (see O Section 4.2]). We 
have Tq > Tq'^° where Tq'^° is the proportion of time of process If spent inside /q. We 
see that T^'^" = Tq'^" where Tq'^" = inf{t > : Yq^° = L,Y/^" =0}. It is not hard to 
prove, via an approximation similar as in [3 Section 4.2], that E^[T(f'^"] = +00. More 
precisely, let 



^ dy i-y 



lo{z)dz 



k{y) Jo dy Jo k{y) 



k{y) 

Then we have E^[T^'^°] = lim v{L;A). Since we can estimate 



A 

dy r 



, , , , lo(z)dz . 



dy 2 Iq 



k{y) 

and we have our Assumption 1, we see that E'^[r,^] = +00. □ 



5 The Large deviation principle 

We are interested in describing the wave front propagation corresponding to the 

solution u{t, (x, k)) of (1.5). To this end we study the quenched large deviation principle 

vt — X^'^^'^Unt) 

for the random variable . Here v > 0, k > and X^'"^'^^{K,t) is the first 

nt ' ^ ^ 

component of the process Yj = {Xt,kt) on F starting from a point {vt,k) G T. Here k 

may be or some other integer > 1 depending on the structure of T. This is in essence 

an adaptation of the arguments of [TH] and [T] . 

Here and below, for notational convenience we will use the symbol X^{Kt) to denote 

the process Xt (which is the first component of the process Yt = (Xf, kt)) starting from 
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a point (x, k) on T with an arbitrary choice of k. The fact that the large deviation 
results for the random variable are independent of the choice of k will 

Kt 

be revealed in the proof of Theorem 5.2. 



Let A € M and we introduce 

g(r,.,A)=E^[e^^'^lT,.<oo] • (5.1) 

Recall that L is the distance between two consecutive vertices Oi and Oj at which 
there is an edge corresponding to a wing. We see that L is a random variable measurable 
with respect to the filtration {J-^}-oo<t<s<oo generated by the shape of D. For each fixed 
shape of D the random variables Tq and Tq^ are well defined and they are measurable 
with respect to the filtration generated by the Wiener process {W^ ,W^). Notice that 
by our Assumption 6 we have cx)>L>L>L>0 where L, L are constants. 

Lemma 5.1. Suppose that A S M is such that 

E(|lnE^[e^^^l^.<J|) <oo ; E (j InE^le^^o'^ J|) < oo . 



Let A G M and c < v. Then almost surely with respect to P the limits 

/,+ (A) ^ lim —L— lnE^[e^^"**l^.*<^] = -^E flnE^[e^^o^l^^^, 

;,-(A) = lim -J_lnE^[e^^- V.<J = -^E flnE^le^^o'^ 

^ ^ ' t^oo{v-c)t ^ T^t<ooi V ^ ^0 ' 

hold. The convergence is uniform with respect to v and c as v and c vary in a set that 
is bounded and {v — c) is bounded away from zero. Moreover, ^^(A) is independent of v 
and c. 

Proof. Let us just work with fj,~^{X). The proof of this fact is essentially the same 
as that of |14l Section 2, Proposition 1] provided we make small modifications. In fact, 
by the strong Markov property of the process on T it is easy to deduce that for 
r < s < t we have 

lng(r, t, A) = lng(r, s, A) + lnq{s, t, A) . 

Let there be located N{n) edges that correspond to the "wings" in the interval x G [0, cn]. 

cn 

We see that lim — — — - = EL holds P-a.s.. On the other hand, we have, by the ergodic 

n-i>oo N[n) 

theorem, that 



In (7(0, cn, A) 
} 

holds P-a.s.. Thus we see that 



hm ^"T;.^"'^^ =E(lnE^[e^^o-l^.,^- 



hm -lng(0,cn,A)= lim ^i^^^l^ = J-E (in E^^le^^o^^ J 



N{n) 
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provided that 

The rest of the argument is the same as in \14:\ Section 2, Proposition 1]. □ 



E(|lnE^[e^^o^^,^J|)<oo. 



We note that by Theorem 4.1 the requirements of Lemma 5.1 always hold for A < 0. 
Let 

= sup{A € M, ^^(A) < oo} . 
Our Theorem 4.1 implies that A^ > 0. 

Lemma 5.2. (Properties of the function ^^(A)) 

(1) /i±(0) = 0; 

(2) //±(A) < /or A < 0; 

(3) fi^{X) — > — oo as A — )■ — oo; 

(4) n^iX) = +00 as X> Xf; 

(5) //^(A) is convex for X € (— oo, A^); 

(6) For X < Xf, /^^(A) is differentiable and 



(/.±)'(A) = E 



In particular, 



(/.±)'(0) = a^ = E[E^[ro±^l^±.<^]] e (0,oo] ; 



(7) {fi^)'{X) is monotonically strictly increasing for X € (— oo, Ac); 

(8) = Oq = +00 and therefore A^ = A~ = 0. 



Proof. The proof of this lemma is the same as in |14l Lemma 2.2, Proposition 
4.2]. The last statement (8) follows from our Lemmas 4.7 and 4.8. □ 



We define 

I^{a) =sup(aA-^^(A)) . 
A<0 

Lemma 5.3. (Properties of the function /^(a)) 

(1) I^{a) >0 forae (0, oo ; 

(2) I^{a) is convex and decreasing in a for a € (0, oo); 

(3) lim /^(a) = +oo and lim /^(a) = 0. 

a— >-0+ a— >-+oo 
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Proof. The proof of this lemma is the same as in |14j . □ 



Theorem 5.1. (Large deviation principle for hitting time) Almost surely with 
respect to P the following estimates hold. Let u,c G M and c < v. For any closed set 
G C (0, oo) we have 

limsup -\nP^ eG] < -{v - c) inf /^ 



t—¥00 



t \ t I aeG \V — C 



limsup -InP^ f — G G ) < -(t; - c) inf F 

t^oc t \ t I a£G \V - C 

and for any open set F C (0, oo) we have 

1 /rpvt 

liminf - InP^ G F 1 > -(v - c) inf / 



t^>oo t \ t J aeG \V — C 

liminf -InP^ (£ F ] > -(v - c) M I" 



t~^oo t \ t I aeG \V — C 



Proof. The proof of this theorem is the same as the proof of \14:\ Theorem 2.3]. 
For the sake of completeness we shall briefly repeat it here. We prove the first and third 
bounds for example. The second and fourth estimates are the same. Let A < 0. Let us 
consider the upper bound first. We have, by Chebyshev's inequality, that 

limsup -InP^ < a 

t— 5>00 t \ t 

<limsup-lnP^ (e^^- >e^"* 
1 

< —Xa + lim sup - In q{ct, vt, A) 

t— s-oo t 



Thus we see that 



= -Aa + {v - c)/i+(A) . 
limsup -InP^ < < - sup(AQ - (t; - c)^+(A)) = -{v - c)I+ ( , 

t^oo t \ t J x<0 \V -Cj 

since A+ = 0. 

We now derive the lower bound. Let u G (0, cxj) and 6 > 0. Let Bs{u) = {u—6, u+5) 
be the (5-ball centered at u. Let A„ < be such that 

/+ (_!^) = sup - ,,+ iX)) = A„^ - /.+ (A„) . 

\v -cJ x<o \ V -c J V -c 

Now we make use of a Cramer's change of measure. Let 
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\ rpvt ^ 



<c»J 



Then we get 
P 



rj-ivt 



e Bs{u) 



One can show m the same way as in [H page 77] and [15], that 

hmmf - InP'^'"'* ^ G S^fn) = . 
Suppose we aheady have (5.2), then we can conchide that we have 



1 / T*'* 

hminf-lnP'^ G Bs(u] 



= {v-c) 
= {v- c)I 
which implies the lower bound. □ 



> -XuU - 5\\u\ + {v - c)^v{\u) 
u 



V — c 
u 

V — c 



-S\Xu\ 



S\Xu\ 



(5.2) 



Theorem 5.2. (Large deviation principle) Almost surely with respect to P the 
following estimates hold. Let v £ U. and n G (0, 1]. For any closed set G C [0, oo) we 
have 



lim sup — In P 



W 



vt - X^Hnt 



and for any open set F C [0, oo) we have 



lim inf — In P 

t->oo Kt 



W 



lim sup — In P 



W 



V 


Kt 


d) we 


have 




-X^HKt) 

^ — - 


^vt - 


Kt 


-00,0] 


we have 


fvt - 


X'"\Kt) ^ 


\ Kt 



€ G ) < - inf cl+ ( - 

CGG V C 



G F I > - inf c/+ I - 



and for any open set F C (— oo, 0] we have 
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li^MllnP"' e f) > - i„f W/- fl 

Proof. We show the first two estimates as an example. The last two estimates are 
the same. We shall make use of the duality 



j^!lll^^!M > = P^ {X^\^t) <Ct- CKt) < P'^ (^S^ < . 

Here = inf{t >0: X'{t) < r}. We have 

= l{Yo /o)^T + T^^ + l(y^. /o)t . 

Here o" is the first time that the process Yt, starting from Yq ^ Iq, arrives at Iq; 
T is the first time that the process Yt arrives at the first branching point K on T with 
x-coordinate > r; r is the first time that the process If, starting from Y^, arrives at 
Yf.^. We note that by our Assumption 6 in probability 1 the distances p(Yo-, {s,0)), 

r 

p{Yt-, {r,0)) is finite. On the other hand, as t — )• oo, by Law of Large Numbers for 
stationary sequences we have — > and — > almost surely. Thus for fixed r < s 

^ rpst rpst 

^ rt ^ rt 

we have 

'J-'St 

^ ^ 1 as i ^ oo . (5.3) 

From here we have 



r 1. ( vt-X-\Kt) ^ \ 1 Tl-c.)t ^ \. 

limsup — mP ^ — - > c < limsup — InP ^ — < k < —cl^ - 

t^oo Kt \ Kt ) t-foo Kt V * / 

which proves the upper bound. 

We now derive the lower bound. We have, for < e < 1, 



P- {'-^^^^ € Bsiu] 



> P^ G ((1 - e)Kt, Kt)) - P^ sup \X^'\s) -{v- Ku)t\ >Kt6] . 

^ ^ \{\-e)Kt<a<Kt j 

The second term in the above formula can be estimated by using space reversal 
invariance and the corresponding large deviation principle, in the same way as |14^ 
proof of Theorem 2.4] and |15t Section 5], provided that we have (5.3). We have 
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lim limsup — InP'^ ( sup \X^^{s) — {v — Ku)t\ > KtS ) = — oo 



e-5>0 t^oo t 



{l — E)Kt<S<Kt 



So then we have, by (5.3) again and Theorem 5.1, 
hmmf-hiP^ ^ G Bsiu) 

t^oo t \ Kt ^ ' 

> Um inf hm inf - In (t^^I^^^, G ((1 - e)^t, Kt)) 
= ( i 

This proves the upper bound. □ 



6 Wave front propagation for reaction diffusion in narrow 
random channels 

After we get the quenched large deviation principle we study the wave front prop- 
agation of the solution u{t, {x, k)) of (1.5) making use of the arguments of [13], [H] and 
[31 Chapter 7]. 

We define non-random constants > and < as the solutions of the equations 

c;/+ = no) , (6.1.1) 

=/'(0) . (6.1.2) 
These solutions exist and are unique due to Lemma 5.3. 

Theorem 6.1. For any closed set F C (— oo,cl) U (c^,cx3) we have 

lim supti(t, (ct, k)) = 

i-i-oo ^^p 

almost surely with respect to P. For any compact set K C {c*_,c*_^) we have 

lim inf u(t, {ct,k)) = 1 

almost surely with respect to P. 



This theorem can be proved in the same way as \14:\ Theorem 1.1, Lemma 4.1, 
Lemma 4.2]. We shall briefly sketch the proof here. We need a sequence of auxiliary 
lemmas. 
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Lemma 6.1. For any closed set F C (— oo,cl) U (c^,+oo) we have 



lim sup u(t, {ct, k)) = 



almost surely. 



Proof. By the KPP condition and (1.5) we have 



{t,{ct,k)) <^lm [exp(/'(0)t)5(X^*(t))] 



We notice that the support of the function g{x) is a compact set U C (—00,00). 
Without loss of generality let us assume that U = Bs = (—6,6) for some 6 > 0. Therefore 
we have 

u{t, {ct, k)) 

<|yexp(/'(0)t)P;rt,fc) H<X^\t)<6) 

n „ N .^W ( ^ Ct-X'^Ht) 5 

= \\g\\ exp(/'(o)t)p;r,,,) + i ^ — ^ " - 7 

We apply Theorem 5.2 with k = 1 and v = c. As t —t- 00 we see that, for 
c > such that /'(O) - c/+ ^-^ < -e < and e > 0; or for c < such that 

/'(O) — |c|/^ ( TT I < — e < and e > we have limsup - lnn(t, (ct, A;)) < — -. This 

\\c\J t-foo t 2 

proves the Lemma. □ 



Lemma 6.2. For any compact set K C (cl, +00) we have 



liminf — In inf n(t, (ct, k)) > — max 

t^oa t ceK ceK 

For any compact set K C (— oo,cl) we have 



c/+ ( ^ 



(6.2.1) 



liminf — In inf uit, (ct, k)) > — max 



(6.2.2) 



Proof. This lemma is proved in the same way as |14l Lemma 4.1], |13l Corollary 
1], provided that we have the estimates (6.3) and (6.4) in the following lemma. □ 

Lemma 6.3. For any y G M and r] > we have 



lim sup P sup \X^{s) — x\ > rjt \ = . 

^->-°° \x\<\v\t \se[o,t] J 



(6.3) 
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Also, for a given M > there exists kq > sufficiently small so that 

limsup sup -InP ( sup - x\ > r]t] < -M , (6.4) 

whenever k < kq. 

Proof. This lemma is proved in the same way as \14:\ Lemma 4.2], by making use 
of Theorem 5.2. We omit the details. □ 

Proof of Theorem 6.1. With the above lemmas at hand the lower bound follows 
from a standard argument as in [13] and |3l Chapter 7, Theorem 3.1]. We omit the 
proof. □ 
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